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, 1 INTRODUCTION 

. Near the Galactic center, stars may get scattered into 
^ ' orbits for which the tide raised by the supermassive black 
, hole, SgrA*, at pericenter is large but not strong enough 
to disrupt the stars. The scattering rate into those orbits 
1 is larger than that of immediate tidal disruptions orbits, 
where the pericenter distances are smaller than the tidal 
Q\ " radius, r v < r t = R*{MbhJM * ) 3 l)Magorrian fc Tremaind 
O ! 1 19991 ; I Alexander fe Livioll200lf). Here M BH = 4 x 10 6 M W i s 
(N ■ the mass of SgrA* jGhez etafl 120081 : iGenzel et al]|20ld h 
and M* and R* are the mass and radius of the star. 
Although stars with r p > rt are not disrupted during 
•"■H ■ their first passage near SgrA*, their tidal heating and 
| bloating could still be substantial after multiple passages 
^ . due to the tidal distortion and the excitation of internal 
■ - - ' oscillation modes. In principle, a sufficiently large number 
of close pass a ges may lead to the d isruption of these stars 



ABSTRACT 

We consider the long-term tidal heating of stars by the supermassive black hole at 
the Galactic center, SgrA*. We show that the accumulated heat deposited by excitation 
of stellar modes over many pericenter passages can lead to a runaway disruption of 
a star at a pericenter distance that is 7-10 times farther than the standard tidal 
disruption radius. The accumulated heating may explain the lack of massive (> IOMq) 
S-stars closer than several tens of AU from SgrA*. 

Key words: black hole physics - galaxies: nuclei - binaries: close - stars: oscillations 




[19881: iNovikov et al 
iDiener et all Il99 



i, 



1 19921: Kosovichev & Novikov 
' Alexander fc Morris! 12003 ; 



Guillochon fc Ramirez- Ruiz! 12012) 1. Various tidal effects 



at r v It r t were considered in the literature, including 
relativistic effects when r t reaches the black hole horizo n 
IlLuminet fc Marckl Il985l : [i vanov fc Chernvako^al I2006T) 
tidal heating of planets by stars l lvanov fc Papaloizoul 
12004a . 2007, 2 0l"lh, and tidal heating in close binary sys- 
tems dPress fc Teukolskvl Il977l; iKochanekl Il992l; iMardling j 

Ivanov fc Papaloizoul 



1995a Ibt lLail 119971: IHo fc Lail 119991: 

2004b1riLai fc Wull2006l; iFuller fc Laill201ll; IWeinberg et al.1 



2012h . 



In this 
distances r f 
ter passage 



paper we consider the heating of stars at 
> 3rt from SgrA*. Since each pericen- 
is associated with a small distortion in the 



shape of stars, one may adopt a linear description for 



the tidal excitation of stella r modes l)Novikov et al 1 ll992l : 
iKosovichev fc Novikovl ll992T ). The associated theory of 
linear mode excitation has been calibrated recently by 
new data on stell a r binaries from th e Kepler satellite 
l|Fuller fc Lail I2OI2I : iBurkart et ail I2012T I. The underlying 
theory was also recently exte nded to describe nonli near cou- 
pling of the excited modes ijWeinberg et al.l l2012h . We use 
the latest results from these studies to calculate the tidal 
excitation and heating of stars in the vicinity of SgrA*. 

Our goal is to find the maximum distance from SgrA* at 
which the accumulated heating due to numerous pericenter 
passages can lead to tidal disruption of stars around SgrA*. 
The accumulated heating would lead to the absence of mas- 
sive stars on eccentric orbits interior to a spherical region 
around SgrA*, whose radius depends on Af* and exceeds the 
standard tidal disruption radius rt- Our predictions could be 
tested by future searches for stars at closer separ ations than 
the known S-stars, whic h have r p > 10 2 AU (|Ghez et al.l 
l200Sl : IGenzel et al.ll2010h . 

S grA* is surrounded by a circumnuclear disk of young 
stars (|Genzel et alj |2010| ). Inside the inner radius of this 
disk, there is the S-cluster of young main s equence B- 
stars (|Ghez et al.ll2003l : lEisenhauer et ahl l2005h ■ with ran- 
dom orbital orientat ions and high orbital eccentricities 
l|Gillessen et al.|[2009r ). All the known S-stars have r p ^> r t , 
but it is possible that the lack of S-stars inside 100 AU is 
caused by the accumulated tidal heating over multiple peri- 
center passages. Our predictions can be tested as new stars, 
such as SO-102 (I Meyer et al.l I2012T ). are being discovered 
and new instruments, s uch as the second- generation VLTI 
instrument GRAVITY (|Bartko et al.ll2009h . are being con- 
structed. 

The outline of the paper is as follows. In Sj2]we describe 
the method we use to calculate the heating due to tidal exci- 
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tation and the response of the stars. In [J3]we show examples 
of these effects in the Galactic center using tw o stellar mod- 
els pr oduced by MESA stellar evolution code (|Paxton et al.l 
l201ll ) and present the results. In [|4] we summarize our main 
conclusions. 



2 HEATING OF STARS BY TIDAL 
EXCITATION OF MODES 

The tidal force from SgrA* can excite internal oscillation 
modes within an orbiting star during its pericenter passages. 
At distances r p > 3r t , the energy gain by tidal excitation 
per pericenter passage is low, but the accumulated energy 
after many passages can heat the star significantly. 



2.1 Mode Excitation and Interference in Multiple 
Pericenter Passages 

To calculate the low energy gain per orbit at r p > 3r t , 
it is appropriate to u s e the linear pe rturbation formalism 
of iPress fc Teukolskvl (119771) (see also iNovikov et al.|[l992l ; 
iKosovichev fc NovikovlTl992l ). We denote the separation of 
the star from SgrA* at time t by r(t). For a single passage, 
the energy of an excited stellar mode can be expressed as, 



A-Eo, 



2tT 



(i) 

where n is the mode order and {I, m} are the two spherical 
harmonic indices. The excited modes have I > 1, — I < m < 
I, and we adopt the convention in which n < for g-modes 
and n > for p-modes. The coefficient K n i m represents the 
coupling to the orbit, 



w lm 

2tt 



dt (^y +iex p{ i [ a;nt + m<i, ( t )]}> ( 2 ) 



where uj n is the mode frequency, &(t) is the true 
anomaly, and W lm = (-1)0+^ ^[^^Q - m )\(l + 

m )!] 1 /2/[ 2 ; ii^li!l!±^!]. The 'tidal overlap integral' Q nl 
represents the coupling of the tidal potential to a given 
mode, 



Q nl = T R 2 dRp(R)lR l ~ 1 [(^ l + + l)e 
Jo 



(3) 



where p(R) is the stellar density profile as a function of ra- 
dius R. £(R) = [£S(i?)e fl + £f i (i?)i?V]Y im (6l,0) is the mode 
eigenfunction, with being its radial component and 
being its the poloidal component. The total energy trans- 
ferred from the orbit to the star in a single passage is 



AE = J2 AEa >> 



(4) 



Next, we consider the evolution of the modes as a result 
of multiple pericenter passages. If the dissipation timescale 
of the modes is longer than the orbital period, the modes re- 
main excited and interfere with newly excited modes during 
subsequent passages. If the orbital period of the star remain 
constant, the mode energy after man y passage s approxi- 
mately equals that in the first passage (Lai 1997). However, 
when the orbital period changes enough (e.g. due to external 



perturbers or the transfer of the orbital energy to mode ex- 
citation energy), the mode a mplitude may increase stochas - 
tically for multiple passages (|lvanov fc Papaloizoul l2004b| y 
Below we show that gravitational scattering on stars and 
compact objects in the Galactic center could naturally lead 
to a drift in the orbital period that allows the amplitude of 
the excited m odes to increase stocha s tically. 

Following llvanov fc Papaloizoul (|2004rJ ). we introduce 
the two-dimensional vectors Xi to characterize the amplitude 
Ai and the phase ipi of the excited modes at the i th passage: 



Ai cos(i/>i), 
Ai sin(V>i)- 



(5) 



Because different stellar modes act independently in the lin- 
ear regime, we focus here on one mode with frequency ui n . 
For the (i + l) th passage, 



U+i = TZ((f>i)[xi + 1 



(6) 



where (pi = u) n P or b,i (with P or b,i being the orbital period 
for the i th passage), e 1 = 1, e 2 = and 1Z is the rotation 
matrix. 

Defining at = ui n AP or b,i where AP or b,i is the change in 
the orbital period in the i th passage, we get 4>i+i — 4>i + a i- 
Given the initial condition xo = (1,0) (without loss of gen- 
erality) and equation ((6]), we examine how the mode am- 
plitude changes a s a function of the number of passages. 
Consistently with llvanov fc Papaloizoul (|2004bl ) . who stud- 
ied the case of ai being a constant, we find that if ai > 0.5, 
the amplitude (|xj|) increases stochastically with i, and the 
threshold ai has a weak dependence on cj>o = uj n Porb,o- For 
Galactic center stars with r p > 3rt around SgrA*, the nec- 
essary change in orbital period per passage cannot be pro- 
vided by the exchange between tidal excitation energy and 
orbital energy. However, it could naturally be accomodated 
by scattering on other stars and compact objects as we show 
next. 

As a prelude to the effect of random scatterings, we 
examined the case when at is drawn from a uniform distri- 
bution between —2a c to 2a c . We characterize the growth in 
the mode amplitude by the power-law index of its evolution 
with the number of passages (using a total of 10 4 passages). 
We aim to find the minimum value of a c above which the 
amplitude increases stochastically. Figure [1] shows that for 
a c > 0.1 the amplitude increases with a power-law index of 
0.5, in analogy to a random walk process in which the energy 
of the mode increases linearly with time. We also examined 
an alternative case with a drawn from a Poisson distribu- 
tion with average (\ati\) = a c and the result was the same. 
In general, we find that the threshold value of a c does not 
show significant dependence on <f>a. 

Next we examine the expected value of a due to grav- 
itational perturbers in the Galactic center. We start by ex- 
pressing a in terms of the fractional change in the orbital 
period assuming the primary excited mode has frequency 
u!„ ~ \f N vm GM* / R% (with a typical value N pm ~ 10), 

a = Ld„AP or b 



- > 3 x 



where e is the orbital eccentricity. Thus, when 
10 -5 the amplitude of the modes increases stochastically. 
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Figure 1. The power-law index of the mode amplitude growth 
with time during multiple passages as a function of the average 
magnitude of a = u>„ AP or (,. When the power-law index is around 
0.5, the amplitude growth resembles a random walk and the en- 
ergy of the mode is growing linearly with the number of passages. 
We find this scaling when a c > 0.1, independent of the value of 
(fio = ^nP O rb,0 (shown by the different lines). 



We calculate the expected AP or {,/P or j, due to grav- 
itational scatterings using th e N-body code BHINT 
l|Lockmann fc Baumgardtl l200ct l to track the orbits of the 
stars and compact objects in the Galactic center. We per- 
formed a convergence test and verified that our numerical 
errors are small and the results are robust. The fractional 
change in the orbital period of a test star depends on the 
semi-major axis a and eccentricity e of its orbit and the dis- 
tribution of perturbers within the S-cluster. We assume an 
outer radius of ~ 0.04 pc (= 1") for the S-cluster, and esti- 
mate the period change for typical S-stars with eccentricities 
in the range of 0.85-0.95. We consider the initial mass func- 
tion (IMF) that matches the mass distribution of S-s tars in- 
side 0.8" (dN/dm oc m- 215±0 ' 3 ) (|Bartko et alJl201Ch ). The 
fractional change of the orbital period is mo st sensitive to 
the massive stars jMurrav-Clav fc Loebll2oTlh . We normal- 
ize the IMF so that it gives ~ 3 S-stars with M ~ 20Mq as 
observed. In the mass range of O.3-25M0, the IMF yields a 
total of 800 stars. 

We also considered the effects of scattering on stellar- 
mass black holes (SBH) and a hypothetical intermediate- 
mass black hole (IMBH). SBHs are more massive than the 
background stars an d therefore are expected to segregate in 
the Galactic center I Morris! 119931 ; IMiralda-Escude fc Gould! 
l200d ; iFreitag et all boOrj ). We normalize the number of 
SBHs (e ach having lOAfo) within 0.0 4 pc to be 1400 
based on IMiralda-Escude fc Gouldl l|2000fl and lFreitag' et al.l 
(2006). An IMBH was hypothesized as an agent for ran- 
domizing the inclinations of stars in the S-cluster, poten- 
tially creating the hyper-velocity sta rs and the stellar disk 
jYu fc Tremaind 120031; ISesana et alj l200d: IVu et all 120071; 
Gualandris fc MerrittJ 120091 ; IPerets fc Gualandris! l20ld ; lYul 
20101 ). To gauge i ts effect on AP or t we assume an IMBH 
mass of 10 3 M Q (YuI |20ld ) with either a = 10~ 3 pc (= 
206AU)) and e = 0.80 or a = 3 x 10" 4 pc and e = 0.26. 



Figure 2. The fractional change in orbital period per pericenter 
passage, AP or i,/P or i,, for stars on orbits with different semi-major 
axis a. The x-axis is in unit of 10 — 3 pc = 206 AU. We include 
800 stars with an initial mass function {dN* / dM*) oc M^T 2 ' 15 
IIBartko et al in the mass range 0.3-25Mq , providing about 

three 20Mq stars. We also consider scattering on a population of 
1,400 stellar-mass black holes (SBH) within 0.04 pc from Sgr A* 
llMiralda-Escude fc Gouldll2000l ; IFreitag et al.ll2006h . or an inter- 
mediate mass black hole (IMBH) with a mass of 1O 3 M0 on two 
possible orbits. 



The scattering due to the SBH and IMBH dominate the 
fractional change in the orbital periods. 

Figure [2] shows the results from the numerical runs of 
the N-body code. We find that the fractional changes of 
the orbital period are higher than the minimum value re- 
quired to increase the mode amplitudes stochastically, im- 
plying that the energy of the excited modes would increase 
linearly with the number of pericenter passages. Because the 
scattering of the orbit is fully random and the change of the 
orbit is typically small (~ 10 _4 -10 -3 ), we neglect the orbital 
evolution. For a random walk, the value of r v is expected to 
change significantly only after 10 6 -10 8 passages, beyond the 
number of passages considered here. 



2.2 Tidal Heating of Stars 

Since the expected fractional change in the orbital period 
per pericenter passage in Figure [2] is higher than 3 x 10 -5 , 
the tidally-excited mode energy is expected to increase lin- 
early with the number of pericenter passages. Cumulatively, 
a significant amount of heat might be deposited inside the 
star during multiple passages. In this section, we consider 
the dissipation of the mode energy and the resultant heat- 
ing of the star. 

Previous studies showed that when the amplitude 
of the excited modes increases over some parametric 
instability threshold, the excited mode begins to transfer 
its energy to lo wer frequency daugh te r modes which dissi- 



pate ra pidly llDziembowskil 19821 ; 

20011; 



Kumar fc Goodmanl 
lArrasetalJ 120031 ; 



,pi< 

19961; IWu fc Goldreich _ 
Weinberg fc Quataerj 20081 ; Weinberg et alj 
set n C rit — Eth/AEo to be the number of pericenter 



20121 ). We 
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passages after which the amplitude of the mode exceeds 
this threshold, where Eth is the threshold energy when 
non-linear coupling occurs. As the dissipation time of the 
excited daughter modes is typically short compared with 
the orbital period in the Galactic center, the thermal energy 
gain in the stellar interior is: 



E, 



= (n p /n cri t)E th = n p AE , 



(8) 



where Et is the thermal energy gained during this process, 
n p is the number of pericenter passages and AEo is the 
energy gain of the excited modes during the first passage. 
When n p 3> n cr u, the thermal energy added to the star is 
independent of the parametric instability threshold. 

The heat generated around a radius R within the star 
at time to will be trapped inside the star for a finite time, 
(t — to) < t c (R), where t c (R) is the characteristic time it 
takes heat to leak out. We estimate t c (R) as the minimum 
between the photon diffusion time, tdiff = / dR{r(R) — 
(R* — R)[dr(R)/dR]}/c, and the turbulent convection time, 
tt = j dR/v c (R), for each spherical shell inside the star. 
Here t(R) is the scattering optical depth and v c (R) is the 
convective velocity. At late times t > t c (R), the heating 
at radius R will saturate and reach a steady state where 
it is balanced by cooling. This sets the upper limit of the 
maximum heat stored at a radius R. 

As the non-linear coupling excites a large number (> 
10 3 ) of daughter modes, most of the energy is redistributed. 
Typically, the daughter modes consist of high order g-modes 
and so the energy is redist ributed mostly in the radiative 
zone. IWeinberg etafl (|2012h investigated modes inside solar- 
type stars and found that most of the energy is transferred to 
the radiative core of the star. For simplicity, we will assume 
that the energy is uniformly distributed per unit mass within 
the radiative zone. 

The energy gained can be expressed as follows, 



E t (R) 



n cri tAEo{R) if t c (R) < P or bn cri t 
ij^-AEoiR) if tc(R) > Por b n crlt 



(9) 



Assuming that energy is evenly deposited throughout the 
entire radiative zone of the star, we find Eo(R) and obtain 
the thermal energy stored at a radius R, Et(R). Typically 
for stars at r p > 3n around SgrA*, n cr it <C t c (R)/P or b, 
and so the total stored heat is independent of n cr u- Finally, 
integrating Et (r) over the interior of the star yields the total 
heating inside the star, Eh- 

As a result of the additi onal source of energy, the star 
expands l|Podsiadlowskill 19961 ). So far, we did not include the 
increase of the size of the star in our calculation. As the stel- 
lar radius increases, the tidal effects become stronger with 
AE oc i?*. A decrease in the mode frequency ~ y/ GM*/R% 
brings ui n closer to the orbital frequency and increases K n i. 
Thus, ignoring the variation in the tidal overlap integral 
(Qni), the tidal excitation becomes stronger as the size of 
the star increases. We examine this pr ocess more closely w ith 
MESA stellar evolution simulations |Paxton et al.ll201lf ) in 
the next section. 

As the star gains energy, its energy gain rate increases 
due to its increasing size. The resulting runaway process 
could lead to the disruption of the star. In order to find the 
minimum heating at saturation (t > t c (R = 0)) that may 
lead to disruption, we express the radius of the star after the 
n th pericenter passage as, 7?»(n) = R*,o(l + e n ), where R t ,o 
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Figure 3. Radius of the star as a function of the number 
of passages after saturation when t > t c (R = 0), assuming 
AR/R*, (n = 0) = 0.1 and 2tt 2 £ ni=2m \Q n i=2\ 2 \K nl=2m \ 2 ~ 
5 X 10~ 3 . We find that a star can be heated significantly after 
~ 10 s — 10 6 passages even if the thermal energy it stores at satu- 
ration is only 10% of its binding energy. This threshold increases 
as r-pjrt increases. 



is the original radius of the star. Assuming AE oc R*(n) 6 
and ignoring the change in entropy within the star, we find 

1/(1 + e„) - 1/(1 + e n+ i) = AE ((1 + en) 6 - 1), (10) 

where AE = AE /(GM* /R*,o)- Figure |3] shows the 
growth of R,(n) as a function of the number of pericenter 
passages starting with the saturation value of eo = 0.1 and 
assuming 2tt 2 £„ !=2m \Qni=2\ 2 \K n i =2m \ 2 ~5x 10~ 3 . Our 
results demonstrate that at r p /rt ~ 5-7 the stored heat can 
approach the binding energy of the star after ~ 10 5-6 peri- 
center passages following saturation, even if the total heat 
gained at saturation is only ~ 10% of the binding energy. 
This threshold increases as r p /rt increases. 

During its lifetime, a massive star can achieve > 10 
pericenter passages at the corresponding distances from 
SgrA*. For example, a 20 Mq star with e = 0.9 and r p ~ 7r t 
around SgrA* has an orbital period of ~ 1.3 years. Thus, 
during its lifetime the star encounters ~ 10 6 pericenter pas- 
sages. The maximum number of pericenter passages is also 
limited by gravitational scatterings on other stars. Accord- 
ing to Figure [2] with stochastic scatterings on SBH or one 
hypothesized IMBH, the maximum number of passages at 
the original pericenter is ~ 10 6 — 10 s . Thus, for r p /r t < 7, 
the star will be significantly heated even if the total heat 
gained at saturation is only ~ 10% of the binding energy. 

Non-linear effects are expected to dominate in the last 
phase of the disruption process. When the star is distorted, 
the energy transfer from the orbit to the modes can be ei- 
ther positive or negative depending on the phases of the 
modes and the orienta tion of the e l lipsoid at the time of 
the pericenter passage. iDiener et alJ ()l995l ) studied this ef- 
fect statistically and found that the probability of a positive 
transfer of energy from the orbit to the star is high. 
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Table 1. Properties of stellar models 



Mass 


Metallicity 


Radius 


Age 


(Mq) 




(Re) 


(yrs) 


1 


Z = Z Q 


l 


4.5 x 10 9 


20 


Z = Z @ 


10 


7 x 10 6 



1.5 



3 1 



US 



x _* xx o r 

mode order (n) 



10 



Figure 4. \Q„ i—2\ 2 \K n 2 m= o| 2 as a function of the mode 
order (n) for stellar modes c omputed with the ADIPLS code 
j Christensen- D aTs gaard 20081) based on the stella r structure 
from the MESA stellar model iPaxton et a.1.1 l20lj]l . The mass 
(20Mp)) and radius (IP-Bp)) of the star resemble those of SO-2 
l|Martins et al. 2008). n < represent the g-modes and ri > 
represents p-modes. 



3 RESULTS 

Based on the formalism presented in we calculated 
the tidal heating of stars in the Galactic center. We con- 
sider two steller masses: IMq (representing low-mass stars) 
and 2OM0 (represen ting high mass stars, similar to SO-2 
(|Martins et al.l l2008)). The other properties of the two stars 
are summarized in Table [1] 

Since the energy gain in each passage depends on 

/ s2i+2 

(R*/r p \ , and because the value of Q n i and K n i m are 
similar for modes with different values of I, the quadrupole 
(/ = 2) modes gain the most energy during the tidal excita- 
tion (whereas I = and I = 1 modes are not excited). Thus, 
we focus on the I — 2 modes. For simplicity, we limit our 
attention to m = modes in non-rotating stars. 

We calculate the overlap integral (Qni) and the 
orbit coupling (K„ im) using the MESA stellar model 
l|Paxton et all l201ll) . The adiabatic normal modes are 
comp uted with the ADIPLS code (|Christensen-Dalsgaardl 
120081 ). For illustration, we show in Figure [4] the values of 
Qn,i=2| 2 |-Kn,i=2,m=o| 2 for a 20Mq star in orbit around Sgr 
A* with a = 5 x 10 -3 pc and e = 0.9. As expected (e.g. 



IPress fc Teukolskvlll977l : iBurkart et al1l2012h . we find that 
lower order modes are excited the most. The energy gain in 
one passage, AEo, can then be found from equation Q. 
To calculate the time it takes for the deposited heat to 




Figure 5. The cooling time (f c ) as a function of radius for the 
two stars in Table [T] 



travel to the surface (t c ) as described in £12.21 we obtain the 
optical depth and the convective veloci ty profile in the inte - 
rior of the stars from the MESA code (jPaxton et al.ll201lh . 
Figure [5] shows the cooling time as a function of radius for 
the two stars. 

The threshold for n o n-line ar coupling has been dis- 
cussed by I Weinberg et ahl (|2012l ) for three mode coupling in 
a solar-mass star. If the daughter modes only couple to one 
other daughter mode, the threshold is E t h ~ 1Q GM 2 /R„ ; 
however, if the daughter modes couple to multiple daugh- 
ters, E th ~ W~ 16 GM 2 /R t . In both cases, (t c /t orb ) > n crlt 
in the interior of the stars for r p > 3r t . Thus, the heating of 
the star is independent of the value of n cr u- Conservatively, 
we calculate the heating using the high energy threshold. 

As the daughter modes consist of high order g-modes, 
energy is redistributed mostly in the radiative zone of the 
star. For simplicity, we assume that the distribution is uni- 
form per unit mass in the radiative zone and integrate equa- 
tion © over the interior of the star. Figure [S] shows the 
heat gained by the stars (Eh) in units of their binding 
energy (Eb) obtained from MESA, at the saturation time 
t — t c (R — 0). The increase in the stellar radius is not in- 
cluded in this calculation. 

Taking account of the runaway increase in the stellar 
radius, the net heat deposited could approach the binding 
energy and hence lead to disruption when the heating at sat- 
uration approaches 10% of the binding energy. We estimate 
that the heating could be substantial at r p ~ 7rt for 20Mq 
stars. 

Next we analyse the heating effect more accurately us- 
ing MESA stellar evolution simulations. We estimate the 
heating rate by AEo/P rb> and assume that the heat is de- 
posited uniformly in the radiative zone. We take account of 
the change in \Q n i=2\ 2 \K n i = 2m\ 2 due to the change of the 
stellar structure through iterations. For our first iteration, 
we assume a constant \Q n i=2\ 2 \K n i = 2m\ 2 and obtain the 
structure of the heated stars with different radii at different 
times. Then we calculate the increase in \Q n i=2\ 2 \Kni=2m\ 2 
as a function of the increase in stellar radius for the heated 
stars. For our second iteration, we simulate the heated stars 
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Figure 6. Maximum amount of heat gained at saturation in units 
of the binding energy of the star as a function of its orbital param- 
eters a and e for a 20Mq star. The solid black line indicates the 
pericenter distance boundary r p = 2.7rt b elow which the linea r 
tidal excitation formalism is not applicable jNovikov et a,l.lll992t) . 
The dashed black line delineates the threshold Eh/Eb ~ 0.1, 
beyond which the star can potentially be disrupted in < 10 6 peri- 
center passages. We find that 20A/q stars are significantly heated 
at r p ~ 7r t . 



10" 



\K„l = 2m | 
2 



as a function of stellar 
2 and continue iterat- 

2 



with a changing |Q„ ;=2 
radius. We calculate \Qni=2f\K n i=2m 

ing until the dependence of \Q n i=2\ 2 \K n i=2m\ 2 on radius con- 
verges. In the examples we consider, convergence is reached 
within two iterations. 

Our convergent results for the 20 Mq star indicate that 
the size of the convective core decreases and the central tem- 
perature stays approximately constant during the heating. 
For the 1 Mq star, the size of the radiative core increases 
and the central temperature drops significantly. Figure [7] 
shows the radius of the heated star as a function of time. 
We compare the results of the two iterations for the 20 Mq 
star at r p /rt = 7 and for the 1 Mq star at r p /r t = 10, 
and find that the disruption time depends only weakly on 
the change in \Q n i=2\ 2 \K n i = 2m\ 2 ■ For other pericenter dis- 
tances we show only the results of the first iteration (as- 
suming \Q n i=2\ 2 \K n i=2m\ 2 =const). Requiring the heating 
timescale to be shorter than the orbital scattering timescale 



and the stellar lifetime, we find that the maximum 



for 



disruption is ~ 7r t for a 20 Mq star and ~ 10rt for a 1 Mq 
star. 

For simp licity, we only consid ered non-rotating stars. As 
discussed bv lFuller fc Lai i|2012r ). the mode frequencies are 
modified for rotating stars by mC n ;fU, where f2* is the ro- 
tation rate of the star and C„j = pR 2 {2( R 'i s +£ S2 ) dR. 
Because Q n i are unchanged by rotation, the dominant 
modes shift to higher order g-modes which have smaller 
values of Q n i. Thus, rotation would lower the excitation 
energies. In addition, the r otatio n may modify the modes 
themselves |Burkart et al.l I2012T ). and further complicate 
the calculation. Treatment of tidal ex citation in misal igned 
spin-orbit system s were discussed by iHo fc Lail i| 19991 ) and 
lLai fc Wul ^OOd ). 
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Figure 7. Stellar radius versus heating time. Top panel: IMq 
star; Bottom panel: 20Mq star. Blue lines indicate the radius 
change by stellar evolution. Requiring the heating time to be 
shorter than the orbital scattering timescale (~ 10 6 yrs) and the 
lifetime of the unheated stars, the maximum r p for which the stel- 
lar radius significantly is ~ 7rt for the 20 Mq star, and ~ 10rt 
for the 1 Mq star. At these limiting cases, the dashed lines show 
results from a second iteration in which \Q n l—2\ 2 \K rl i—2m\ 2 ls 
updated as the stellar radius increases. 



Finally, we discuss the observational signature of a 
tidally heated star. Using the MESA simulation, we plot 
the Hertzsprung-Russell (HR) diagram of the heated stars 
in Figure [S] Because our calculation is not appropriate in 
the non-linear regime when the tidal radius of the heated 
star approaches ~ (r p /2.7), we stop the calculation when 
r p ~ 2.7r ttheated , where r t , heated is the tidal radius of the 
heated star. We find that a 1 Mq star at r p ~ 2.7rt, heated 
acquires a luminosity L that is ~ 6 times higher than if 
it were on the main sequence and an effective temperature 
T e ff that is ~ 20% lower than the main sequence star. A 
20M© star at r p ~ 2.7rj. heated acquires a luminosity that 
is ~ 2 times higher and an effective temperature that is 
~ 33% lower than on the main sequence. Photometrically, 
the heated stars could be confused with giant stars that 
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Figure 8. HR diagram of heated stars with masses of IMq (left 
panel) and 20Mq (right panel). Blue lines indicate the evolution 
track of giant stars with the same masses as they evolve off the 
main sequence. The HR diagrams of the heated stars stop at the 
point when the tidal radii of the heated stars approach (r p /2.7), 
at which point the linear tidal excitation approach breaks down. 



evolved off the main sequence (illustrated by the blue lines 
in the plot). 



4 CONCLUSIONS 

We considered the tidal excitation of oscillation modes in 
stars orbiting SgrA*. When the dissipation timescale of the 
modes is longer than the orbital period, the modes excited in 
each passage interfere. Due to the gravitational scatterings 
on nearby stars or stellar-mass black holes, the orbital period 
of the excited star changes stochastically and the energy of 
the excited modes increases approximately linearly with the 
number of pericenter passages. As non-linear coupling of the 
stellar modes dissipate the kinetic energy of the modes, the 
excited star is heated. Once the deposited heat is significant, 
the star bloats and its tidal heating accelerates, until non- 
linearities lead to the final mass loss and possible disruption 
of the star. 

We calculated the thermal energy gain by a star as a 
function of the semi-major axis and eccentricity of its orbit 
around Sgr A*. We have found that the maximum pericenter 
distance where the heat gained by the star approaches its 
binding energy is r p ~ 10r t (~ 7.5 AU) for a 1 M@ star 
and r p ~ 7r t (~ 20 AU) for a 20 Mq star. The accumulated 
heating may explain the lack of massive (> I OMq) S-stars 
close r than several tens of AU from SgrA* (Genz el et al.l 
120101 ). 

The heating process may be most effective for the 
highest-mass stars (> lOOA/ ), where radiation pressure 
nearly balances gravity and red uces the binding energy con - 
siderably relative to GMt/R, (|Shapiro fc Teukolskvlll98rj) . 
This makes these stars more vulnerable to disruption 
through heating. However, the heating is not important 
for giant stars evolved off the main sequence, because for 
fp/rt ~ 10 the orbital period of a giant star is too long 



to allow sufficient number of pericenter passages during the 
star's lifetime. 

The expected radius of the cavity pr oduced by tidal dis- 
rupti on of stars depends on stellar mass ^Alexander fc Livid 
2001). Since gravitational scatterings on other objects could 
change the value of r p on a timescale much shorter than the 
lifetime of a low mass star but similar to the lifetime of the 
high mass star (~ 20Mq), the net number of pericenter pas- 
sages for a given r v is similar in the two cases. Of course, 
the tidal distance of a high mass star is larger than that of 
a low mass star, and so a lower mass star may approach Sgr 
A* at a closer distance (having a shorter orbital time and 
more pericenter passages) before being tidally disrupted. 

The removal of tidally-heated stars makes it more diffi- 
cult to test the no hair theorem of general relativity based 
on stellar orbits, as the precession produced by the quadru- 
ple moment of SgrA* decreases with increasing distance. For 
example, the precession rate due to the quadruple moment 
of SgrA* is only ~ 0.1/ias/yr for a 20Mq star with r v — 7rt, 
and is ~ 1/ias/yr for a IMq star with r„ = 1 0rt, assuming a 
normalized spin of 0.7 for SgrA* (I Willi 120081 ). Gravitational 
deflections by other stars or compact objects contaminate 
the precession signal and require t he monitored stars to be 
within ~2x 10" 4 pc from SgrA* (|Merritt et alJl201Ch . We 
find that only low mass stars (which cannot be detected 
at present) would be viable targets for testing the no hair 
theorem around SgrA*. 

As new stars, such as SO102 (jMever et al.ll2012r i. are 
being discovered in the Galactic center, our predictions for 
the tidal cavity radius as a function of stellar mass may 
be tested. In particular, the second-generation VLTI instru- 
ment GRAVITY will be able to resolv e faint stars with a 
K-band magnitude m K = 18 (~ 3M ) (|Bartko et al.ll2009f ) 
and test our predictions in the coming years. 
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